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ABSTRACT 

Recent observations have revealed that damped Lya clouds (DLAs) host star formation activ- 
ity. In order to examine if such star formation activity can be triggered by ionization fronts, 
we perform high-resolution hydrodynamics and radiative transfer simulations of the effect of 
radiative feedback from propagating ionization fronts on high-density clumps. We examine 
two sources of ultraviolet (UV) radiation field to which high-redshift (z ~ 3) galaxies could 
be exposed: one corresponding to the UV radiation originating from stars within the DLA, 
itself, and the other corresponding to the UV background radiation. We find that, for larger 
clouds, the propagating I-fronts created by local stellar sources can trigger cooling instability 
and collapse of significant part, up to 85%, of the cloud, creating conditions for star forma- 
tion in a timescale of a few Myr. The passage of the I-front also triggers collapse of smaller 
clumps (with radii below ~ 4 pc), but in these cases the resulting cold and dense gas does 
not reach conditions conducive to star formation. Assuming that 85% of the gas initially in 
the clump is converted into stars, we obtain a star formation rate of ~ 0.25 Mq yr _1 kpc~ 2 . 
This is somewhat higher than the value derived from recent observations. On the other hand, 
the background UV radiation which has harder spectrum fails to trigger cooling and collapse. 
Instead, the hard photons which have long mean-free-path heat the dense clumps, which as 
a result expand and essentially dissolve in the ambient medium. Therefore, the star forma- 
tion activity in DLAs is strongly regulated by the radiative feedback, both from the external 
UV background and internal stellar sources and we predict quiescent evolution of DLAs (not 
starburst-like evolution). 

Key words: hydrodynamics — radiative transfer — methods: numerical — ISM: molecules 
— galaxies: evolution — quasar absorption lines 



1 INTRODUCTION 

Damped Lya clouds are quasar (QSO) absorption line systems 
whose neutral hydrog en column dens ity is larger than ~ 2 x 10 20 
cm~ 2 (e.g. Prochaska & Wolfe 2002). This class of objects dom- 
inates the neutral hydrogen content aroun d z ~ 3 a nd possibly 
also at higher and lower redshifts fe.g. lPeroux et all2003l) . Because 
QSOs are generally luminous, DLAs provide us with unique oppor- 
tunities to trace high-redshift (high-z) galaxy evolution. In particu- 
lar, by identifying absorption lines of various species, the physical 
condition of the interstellar medium (ISM) of DLAs has been de- 
duced. 

The detection of metal lines in DLAs demonstrates the exis- 
tence of heavy elements in DLAs fe.g. lLu et alJl99rilLedoux et alJ 
2003). The evolution of metal abundances in DLAs can trace the 
chemical enrichment history of present day galaxies. Based on this 
and other clues, DLAs have been suggested to be the progenitors 
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of nearby galaxies; the similar values of the baryonic mass den- 
sity in DLAs around redshift z ~ 2 and the stell ar mass den- 
sity at z ~ has further support ed this idea lLanzetta et all 19951 : 
IStorrie-Lombardi & W olfe 2000). By adoptin g a recently favoured 
Lambda cold dark matter cosmology, however, Peroux et aljfcoojl) 
showed that the comoving density of H I gas at z ~ 2 is smaller 
than the comoving stellar mass density at z ~ 0. Yet they empha- 
size that a large fraction of H I gas is contained in DLAs at z ~ 2-3. 

The existence of heavy elements in DLAs indicates that DLAs 
have experienced star formation. There is also direct evidence of 
ongoing star formation activity from the UV interstellar radia- 
tion field larger than the ultraviolet b ackground (UVB) intensity 
jLedoux et all2002tlWolfe et alj2003h . lHirashita etai]<2003l) have 
theoretically shown that there are small-scale clumps whose molec- 
ular fraction is higher than 10 -3 in DLAs. The typical size and den- 
sity of such clumps are ~ 5 pc and 30 cm" , respectively. How- 
ever, it is not clear if those clouds finally collapse and form stars 
in UV radiation field. It is however crucial to address this question, 
as most these clumps are likely to host the most active sites of star 



2 I. T. Iliev, et al. 



formation in DLAs. iHowk et"aT] 120051) have shown the existence 
of cold gas in a high- 2 DLA, arguing that such cold gas may work 
as the fuel for star formation. 

The nature of DLAs is at present uncertain. They have been 
thoug ht to be the p rogenitors of nearby large galaxies. For example, 



Wolfe et al. 1 198 (j) propos e that DLAs are disk galaxies ( see also 



Prochaska & Wolfe 1998). More recently, Haehnelt 



des (s < 
et all 



and Ledoux et al. 1 1998) showed that the velocity structure can be 
explained by infalling subgalactic clumps in collapsing dark matter 
halo. At z ~ 0, spiral gala xies may c ontribute predominantly to 
the damped Lya absorption l Rao & Briggs 1993), but a t high red- 
shift, various classes of objects can be observed as DLAs ICen et ail 
2003). In any scenario, small clu mps can form as a resu lt of non- 
linear hydrodynamical processes l Wada & Nor marl200ll) . and it is 
important to know if such clumps finally lead to star formation ac- 
tivity. 

In a previous work l Hirashita_et_alJ|2003j), we h ave used the 
2-dimensional (2-D) code of Iwada^^Normanl 1200 ll) suitable for 
disk-like geometry, in order to follow the evolution of entire galac- 
tic disk. They stress the clumpy structure of gas, whose typical 
size is roughly ~ 5 pc. In this paper, we build on this by examining 
the fate of the resulting clumps in external UV field. For this prob- 
lem, we use a 2-D axially-symmetric radiative-hydrodynamical 
Adap t ive M e sh Refinement (AMR) code CORAL iRaga et alJ 
I l995t Il997t lMeIlemaetalJll998t IShapiro. Iliev & Rag A 1200? 
Illiev. Shapiro. & Raga 1 120051) . which is appropriate for treating 
each clump. 

The structure of this work is as follows. We describe our nu- 
merical methods and input parameters in § |2| In section § [3] we 
present and discuss our results and finally, in §|4]we give our sum- 
mary and conclusions. Throughout the paper we use CGS units, 
unless otherwise noted. 



2 METHOD 



2.1 Numerical methods 



We use a 2-D axially-symmetric radiative-hydrodynamical Adap- 
tive M esh R e finement (AMR) Eule r ian code CORAL jRagaetalJ 

E ll997t lMellemaetalJll998t IShapiro. Iliev & Rag A 120041 
Shap iro. & Raga 1 120051) . The hydrodynamic equations are 
solved by a van Leer flux- splitting method Ivan LeeJ fl982). im- 
proved to second orde r accu r acy by use of linear gradients within 
cells as described in lArthud Il99lh . The microphysical processes 
are coupled to the hydrodynamical evolution by operator-splitting, 
i.e. solved at each timestep and results included in the energy 
equation. These microphysical processes include photo- and colli- 
sional ionizations, recombinations (both radiative and dielectronic), 
charge exchange and transfer of ionizing radiation. The code fol- 
lows the non-equilibrium chemical evolution of the ionic species of 
H, He, C II-VI, N I-VI, O I-VI, Ne I-VI, and S II-VI (C I and S I are 
assumed fully-ionized). Self-gravity is currently not included. The 
refinement (and de-refinement) of a given region in space is decided 
based on the gradients of all code variables, when the gradient is 
larger than a pre-defined value, the region is refined and vice- versa. 
We used a total of 8 levels of refinement (separated by factors of 2 
in cell size). More details on I-front tracking implementation and 
detaile d tests of the method can be found in Sha piro. Iliev & Ragal 
(2004). 



2.2 H2 formation and destruction 

For the metallicity range typical of DLAs, we can assume equilib- 
rium between H2 formation and destruction, because the timescale 
of H2 formation and destruction is well below the dynamical 
timescale (Hirashita et al. 2003). We adopt the formation rate of H2 
per unit volume and time, -Rd ust, bv lHollenbach & McKed 1 19791) 
(see also Hirashita et al. 20Q3) 1 : 
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where a is the radius of a grain (assumed to be spherical with a 
radius of 0.1 /im in this paper), T> is the dust-to-gas mass ratio, S is 
the grain material density (assumed to be 2 g cm -3 in this paper), 
and Sd(T) is the sticking coefficient of hydrogen atoms onto dust . 
The sticking c oefficient is given by iHollenbac h & McKedl 19791 
Omu kail200a) 

S d (T) = [1 + 0.04(T + T d ) - 5 + 2 x 10 _3 T + 8 x lO^T 2 ] -1 
x [1 + exp(7.5 x 10 2 (l/75 - 1/T d ))] _1 , (2) 

where Td is the dust temperature. We adopt Td = 20 K, a typ- 
ical dust temperature in equilib rium wit h the Galactic interstellar 
radiation field (ISRF) I Shib ai et alll999l) . However, since the reac- 
tion rate is insensitive to the dust temperature as long as Td ^ 70 
K, the following results are not affected by the dust temperature. 
In fact, some D LAs show ISRF stren gths similar to the G alac- 
tic va lue <Ge & BechtolJl997tlLedoux et all2002l:|p"etitiean et alJ 
2002). The H2 formation rate per unit volume is estimated by 
Ri = Rdustnnn(l — /h 2 )> where n — 1.08nH is the total, hy- 
drogen and helium number density, uh is the hydrogen number 
density and /h 2 = 2nH 2 /^H (jih 2 is the H2 number density) is 
the H2 molecule fraction. 

T he photodissociation rate J?di ss is estimated by jAbel et al 
1997) 



(3) 



^diss = (47r) 1.1 x 10 s Jlw Shield s 1 ■ 



where Jlw (erg s _1 cm -2 Hz -1 sir -1 ) is the l/Kintensity at hu — 
12.87 eV averaged for all the solid angle, Shield is the correction 
factor of the reaction rate for H2 self-shielding and dust extinc- 
tion. The photodissociation rate per unit volume is then estimated 
by R2 = -Rdiss?iH2, where riH 2 is the number densit y of H2. We 
adopt the correction for the H2 self-shielding by Drai ne & Bertoldil 
( 1996)(see also Hirashita & Ferrara 2002). Then, we estimate 

Sshield aS 
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where Nh 2 is the column density of H2. The last term, which is due 
to absorption by dust is in fact of no r elevance for the parameters 
we assume here iHirashita & Ferraral2005l) . 

The metallicity level of DLAs, ~ 0.1 Z Q IPettini et alll994l) . 
which is the value we assume here, implies that the dust-to-gas ratio 
of DLAs is typically 10% of the Galactic (Milky Way) value. In 



1 In iHirashita et all 120031) there is a typographic error. (The results are 
correctly calc ulated.) The expression in this paper is correct and R\ in 
IHirashita et alJ l2003t) is equal to i?dust nrl H(l — /h 2 )> where /jj 2 is the 
molecular fraction. 
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Figure 1. Distributions of clumps in disk in (a)(left) mean density and (b)(right) size. 
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Figure 2. Evolution of (top to bottom panels): the position of the I-front along x-axis, xj, the neutral mass fraction remaining in the clump, Mi /M- in n, the 
velocity of the I-front vj, and the ratio of the cold, neutral mass to the current (averaged) Jeans mass, Mj, for Cases 1 (left), 2 (middle) and 3 (right). 



this paper, we assume the Galactic dust-to-gas ratio to be 0.01. We 
define the normalized dust-to-gas ratio k as 

V = kV q , (5) 

where T>q is the Galactic dust-to-gas ratio (i.e., Dq — 0.01), thus 
we assume k = 0.1 and hence V — 10 -3 . 

The typical Galactic ISRF intensity has been estimated to 
be cvu v = 1.2 x 10~ 3 erg cm -3 at the wavelength of 1000 A 
(i.e., v = 3.0 x 10 15 Hz ), where u v is the energy density of 
photon per unit frequency (Habina 1968). Approximating the en- 
ergy density of photons at 1000 A with that at the Lyman-Werner 
Band, we obtain Jlw at the solar vicinity, Jlwo, is approximately 
J LW0 ^ cu„/4tt = 3.2 x 10~ 20 erg cm" 2 s" 1 Hz" 1 sr" 1 . The 
intensity normalized by the Galactic ISRF, x< is introduced by 

Jlw = xJlwq ■ (6) 
Using equations {3J and |6|, we obtain 

R diss = 4.4 x lO-^x&hicid s" 1 . (7) 



When x — 1» me typical Galactic photodissociation rate (2-5 x 
10" 11 s~ x ) derived bv ljuriJl 1974b is recovered. 

The UV background intensity at the Lyman limit is taken to 
be J21 = 0.3-1 around z ~ 3, where J21 is in units of 10 -21 
erg cm -2 s _1 Hz -1 sr" 1 (Giallongo et al. 1996; Cooke, Espey, & 
Carswell 1997; Scott et al. 2000; Bianchi, Cristiani, & Kim 2001). 
For the UVB, we assume that J21 = 0.6 and that the spectrum of 
ionizing radiation is described with a power law with slope of 1.55. 
We also take into account a "boost" of 10 1 ' 5 of the Lyman-Werner 
photons flux with respect to the ionizing flux at the threshold, ac- 
cording to the results of Haardt & Madau (1996) at z — 2 — 4 (i.e. 
we assume x = 0.6). 

The UV radiation field originating from stars within DLAs 
may be larger than the UVB (Wolfe et al. 2003; Hirashita & Fer- 
rara 2005). In particular, the UVB intensity is likely to be lower at 
2 £ 1 than at z ~ 3 (e.g. Scott et al. 2002). Therefore, we also 
examine the case where the UV field is dominated by internal stars. 
Hirashita & Ferrara (2005) have derived x = 3 — 30 by using obser- 
vational data of H2-detected DLAs. We adopt x — 15 in this paper. 
Although it is known that ionizing photons are strongly absorbed 
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Table 1. Simulation parameters 



Quantity 


Case 1 


Case 2 


Case 3 


spectrum 


BB 5e4 


BB 5e4 


PL, a = 


^clump (pc) 


5 


3 


5 


"clump (cm -3 ) 


30 


30 


30 


"out (cm -3 ) 


0.3 


0.3 


0.3 


-^clump (K) 


100 


100 


100 


T out (K) 


10 4 


10 4 


10 4 


J21 


6 


6 


0.6 


Mcloud (%) 


474 


102 


474 


/H 2 ,init 


lO" 8 


lO" 8 


lO" 8 



by the neutral hydrogen within DLAs, there is little quantitative 
constraint on the intensity of ionizing photons. 

However, it would be reasonable to assume that 90%-99% 
(Ciardi, Bianchi, & Ferrara 2002) of the ionizing photons are ab- 
sorbed relative to the dissociating photons in passing through the 
interstellar medium, the "boost" could be 10-100. Assuming the 
boost of 50, we obtain J21 = 6 for the internal source. The spec- 
tral shape of the 50,000 K black-body is assumed for the internal 
source. 

In summary, in our simulations we fix the value of the external 
ionizing flux entering the computational box, F, to be equivalent to 
a mean isotropic background J of either J21 = 6, with 50,000 K 
black-body spectrum for the cases where the ionizing sources and 
internal, nearby hot stars, or of J21 = 0.6, with power-law spec- 
trum with slope of 1.55, for the mean intergalactic background (or 
possibly internal/nearby QSO sources). The level and spectrum of 
the ionizing flux F(i,j) reaching each computational cell (i, j) is 
determined self-consistently by the code based on the optical depth 
from the source to the cell. Finally, the intensity in the Lyman- 
Werner bands (normalized to the Galactic value), x> is assumed to 
be 0.6 for the UVB and 15 for the internal sources. 

As mentioned at the beginning of this subsection, we can as- 
sume that the formation and destruction of H2 are in equilibrium. 
Therefore, the following equation holds: 



Rdustnnn(l — /h 2 ) = Rdi ss nH 2 ■ 

Using equation {5}, the molecular fraction is obtained as: 

2nH 2 2i?dustn 
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3 RESULTS 

We have performed three high-resolution simulations with varying 
cloud sizes and external ionizing flux level and spectrum, summa- 
rized in Table Q Cloud sizes and densities were chosen similar to 
the ones found on small sc ales in the larger-scale simulations of a 
galactic disk performed by Hira shita et alJ |2003). In FigureQwe 
show the density and size distribution functions of the clouds from 
the simulations in iHirashita et ail {2003 ), where each cloud is iden- 
tified as a 2-D region with local gas density exceeding 10 cm -3 
and the cloud radius is defined as (A/n) 1 ^ 2 (A is the area of the 
cloud). 

In all cases of our simulations listed in Table 1, both the 
clump and the external medium start neutral, and the clumps are 



in a pressure-supported equilibrium. Our fiducial case (Case 1) is a 
spherical cloud of radius r — 5 pc, temperature T = 100 K and gas 
number density of n — 30cm~ 3 (corresponding to a total clump 
mass of M — 474 Mq), and surrounding medium with tempera- 
ture T — 10 4 K and gas number density of n = 0.3 cm . The 
external flux has 50,000 K black-body spectrum and J21 = 6. Our 
Case 2 simulation has all the same initial conditions as Case 1 , ex- 
cept the clump is smaller, with initial radius of 3 pc (corresponding 
to a total clump mass of M — 102 Mq). Finally, we considered a 
case with power-law external ionizing flux, / oc with a slope 
of a = 1.55 and J21 = 0.6, and all other initial conditions the 
same as in Case 1 . The Lyman- Werner flux is fixed in all cases by 
Jlw = 474 J21, as described in the previous section. 

The Inverse Stromgren length, is, i.e. the size along a line-of- 
sight of the inverse static Stromgren volume from the source side 
into the clump gas along the axis of symmetry, is given by 

= ~WT> < 10 > 

whe re off is the Case B recomb ination coefficient for hydrogen 
[see IShapiro, Iliev & Raeal l2004t) for details]. Numerically it is 
equal to £ s = (0.005,0.005,0.0005) pc in our Cases 1, 2 and 
3, respectively, much smaller than the corresponding cloud diame- 
ter. Therefore the I-front would always convert to slow, D-type and 
become trapped by the cloud. 

Cases 1 and 2 were run with 1024 x 2048 finest grid 
resolution, while Case 3 with 512 x 1024, since simulations 
which adopt harder ionizin g spectra converge more readily 
llliev. Shapi ro. & Ragall2005l) . and thus do not require such high 
resolution. We have also tested our results for convergence by run- 
ning lower-resolution simulation of Case 1, at 512 x 1024 finest- 
grid resolution. Our results show almost no dependence on reso- 
lution; thus our conclusions remain unchanged and unaffected by 
resolution concerns. 

We show the evolution of the I-front position, velocity, and the 
neutral mass fraction remaining in the clump for Cases 1, 2 and 3 in 
Figure [2] We note that there are significant differences in the evo- 
lution of these quantities in the three cases. In Case 1 about 15% of 
the cloud mass becomes evaporated within the first 3 Myr after the 
arrival of the I-front, and no mass is lost afterward, with the cloud 
retaining the rest, about 85%, of its initial mass. The I-front initially 
propagates into the cloud with vi ~ 2 km s ~~ 1 . After that, its veloc- 
ity continually decreases until the I-front stalls at t ~ 3 Myr. The 
cloud itself is also moving during that time at ~ 0.1 — 0.2km s _1 , 
a manifestation of the well-known rocket effect, whereby the pho- 
toevaporating gas pushes the cloud away from the source due to 
conservation of momentum. In Case 2 the evolution is initially sim- 
ilar to the one observed in Case 1, with the cloud losing ~ 25% of 
its initial mass in the first 3 Myr due to photoevaporation with a 
subsequent decrease in the mass loss rate. However, in this case the 
outflow still continues afterward, albeit at a low level (the sharp de- 
crease of the neutral mass fraction after 3.5 Myr is largely due to 
the cloud partially leaving the box, rather than to mass outflow, as 
would be seen below). The clump is again pushed away from the 
source due to rocket effect, at somewhat larger speed (~ 2 kms -1 ). 
Finally, in Case 3 the cloud stays neutral for large part of its evo- 
lution, expanding into the surrounding gas until it eventually be- 
comes mostly ionized between 3 Myr and 7 Myr after the arrival of 
the ionization front and disperses. 

These cloud collapse/dispersion timescales due to radiative 
feedback, which are of order 3-7 Myr, are similar to the forma- 
tion/survival timescales obtained from the gas-dynamic Simula- 
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Figure 3. Case 1 simulation results at (a) t = 2 Myrs, (b) t = 3 Myrs and (c) t = 4 Myrs. From top to bottom: (1) grayscale image of the decimal logarithm 
of the atomic density in (r, a;)— plane of cylindrical coordinates; (2) flow velocities; arrows are plotted with length proportional to gas velocity. An arrow of 
length equal to the spacing between arrows has velocity 25 kms" 1 (i.e. arrows are touching along x- and y-axes if v = 25 km s — 1 ); the minimum velocities 
plotted are 3kms _1 . Solid line indicates the boundary of the gas which was originally inside the cloud, dashed line indicates the current position of the 
ionization front (50% ionization of H); (3) isocontours of H2 fraction, /jj 2 , logarithmically-spaced (from 1 down by factors of 10); cuts along the r = axis 
of: (4) molecular hydrogen column density, Nh 2 , (cm -2 ) and (5) temperature (K). 



tions in lHirashita et alJ 120031) . which can be roughly estimated as 

Survive ~ r c i oud /«flow ~ 5pc/lkms _1 = 5Myr, where « flow 
is a typical gas flow velocity observed in the simulations. In the 
following, we discuss in more detail each of the Cases. 

3.1 Case 1 

Figure [5] shows snapshots of the gas flow for Case 1 at t = 2, 3 
and 4.5 Myr. When the I-front reaches the clump (at time a frac- 
tion of a Myr) it slows down considerably, converts to D-type and 
gets trapped. Initially, as the I-front is entering the clump, some 
Raleigh-Taylor instabilities develop at the boundary between the 



hot gas outside and the cold gas inside. The I-front is clearly pre- 
ceded by a (oblique) shock (at the x-axis it is at x ~ 6 x 10 19 cm 
at time t — 2 Myr), as expected for a D-type I-front. A fraction of 
the cloud material is photoheated and evaporated, outflowing in a 
supersonic wind with speeds ~ 25kms _1 . The rest of the clump 
gas is strongly compressed by the propagating shocks. The dense, 
shocked neutral gas has sufficient column density to partially shield 
the unshocked part of the cloud and the hydrogen is already al- 
most fully-molecular /h 2 ~ 1 there, while the H2 column density 
through the cloud, -/Vh 2 > exceeds 10 19 cm -2 , and the molecular 
hydrogen self-shields. By t — 2.5 Myr the right-propagating shock 
which precedes the I-front crosses the whole clump. At t = 3 Myr 
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Figure 4. Same as in Fig.0 but for Case 2 at times (a) t = 1 Myr, (b) t =2 Myr, and (c) t =4 Myr. 



(Figure [5] (b)) the clump is compressed to less than a third of its 
original radius, and corresponding gas densities n £ 10 3 cm -3 . By 
this time the neutral portion of the cloud is fully self-shielded and 
the hydrogen is completely molecular (/h 2 ~ 1) there, while the 
column density of H2, Nu 2 . exceeds 10 22 cm -2 . By t — 4.5 Myr 
the photoevaporation wind almost disappears, with the neutral mass 
of the cloud remaining stable. The clump is again close to uniform 
at this time, but compressed to a fraction of its original volume and 
gas number density of n ~ 10 3 cm -3 . The clump remains self- 
shielded, contains about 85% of its original gas, and the hydrogen 
inside is largely in molecular form. The Jeans mass is given by 

/ U„ \ 3/2 T 3/2 

where ks is the Boltzmann constant and G is the gravity constant 
and p is the gas mass density. In Figure |2| we plot the ratio of the 



mass of the cold, neutral gas remaining to the mean Jeans mass at 
that time for all three simulations. The average Jeans mass is cal- 
culated based on the mass-weighted temperature and the mean gas 
density. For the compressed cloud (t ~ 4 Myr) we have approx- 
imately T ?s 100 K and n w 500 cm" 3 , and Mj w 377 M©, 
while the mass of the neutral gas in the cloud at this point is 
about 4OOA/0 , hence the cloud should become Jeans-unstable (see 
Fig- 13- Later on (t ~ 4.5 Myr) the cloud expands and the den- 
sity drops slightly, but the temperature also drops due to adiabatic 
cooling to T ~ 60 K, and the local Jeans mass decreases a bit 
more. The physical conditions reached suggest the cloud becomes 
gravitationally-unstable, cold and molecular, which should lead to 
increased star formation in this case, triggered by the original I- 
front. 
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3.2 Case 2 

Results from our Case 2 at i = 1,2 and 4 Myr are shown in Fig- 
ure01 In this case again the I-front is trapped in the cloud and con- 
verted to a D-type, preceded by a shock, and expanding a super- 
sonic wind into the surrounding gas, as is clearly seen in Figure[4] 
At t ~ 1 Myr, some molecular hydrogen is created in the pre-shock 
part of the clump, but its fraction is very small, less than 10~ 4 . Nh_ 2 
reaches ~ 10 17 ' 5 cm -2 , and the cloud self-shields. 

The shocks again compress the cloud, similarly to Case 1, but 
the maximum gas number density reached in the course of the evo- 
lution is somewhat lower than in Case 1, never exceeding 10 3 cm -3 
(Figure|4j. At t ~ 3 Myr the clump reaches maximum density and 
the molecular fraction is /h 2 ~ 1, while Nu 2 reaches maximum 
of ~ 10 20 ' 5 cm~ 2 , becoming completely self-shielded. 

Afterward, the gas outflow due to photoevaporation contin- 
ues, albeit at a somewhat diminished level and the remaining cloud 
is accelerated. As the cloud is pushed back it expands and elon- 
gates, and its density decreases ton ~ 10 2 cm~ 3 . By t — 4 Myr 
much of the molecular hydrogen is exposed to the radiation again 
and mostly destroyed, with /h 2 decreasing to less than 10~ 2 every- 
where, while 7Vh 2 decreases to ~ 10 18 cm -2 , however the cloud 
is still self-shielded. 

At the point of maximum compression its density and temper- 
ature are n ~ 10 3 cm -3 and T ~ 100 K, respectively, and the 
local Jeans mass exceeds the mass of neutral gas remaining in the 
cloud. After the cloud expands again (t ~ 4 Myr) its density drops 
to n ~ 10 2 cm -3 and its temperature drops to T ~ 60 K, in which 
case the Jeans mass becomes Mj ~ 600 Mq, while the neutral 
gas remaining in the cloud is M ~ 70 Mq, i.e. it never becomes 
gravitationally-unstable (see Fig.[2j- While the cloud is still self- 
shielded at that time, the hydrogen molecule fraction is low: and 
conditions are not conducive to star-formation in this case. 

3.3 Case 3 

The results for Case 3 are illustrated in Figure [5] The evolution 
in this case is quite different from the previous two cases. The I- 



front is initially trapped outside the clump itself, but the hard pho- 
tons which are present in this case reach the cloud, since they have 
much longer mean free path and start heating it. The cloud becomes 
heated to T ~ 10 3 — 10 4 K and never gets strongly compressed, 
but instead expands almost uniformly in all directions due to the 
higher pressure inside it. It eventually evaporates completely and 
never forms significant molecular hydrogen fraction. 



3.4 Comparison with observations 

Recent observations have revealed that DLAs host star formation 
activity (e.g., Wolfe et al. 2003). We have found that propagating 
I-fronts produced by internal stellar sources can trigger cooling 
instability (i.e. a "compression-cooling-further compression" cy- 
cle) and collapse of average clumps (with r ^ 5 pc and n ~ 30 
cm _3 )in DLAs. The hydrogen in the collapsed clumps is largely 
in molecular form and the temperature stays below 100 K. The 
physical conditions are, therefore, favourable for star formation. 
This star formation is triggered by the gas compression due to the 
passage of the I-front. Hence, the UV produced by internal stellar 
sources works as a "positive feedback" for the star formation in 
DLAs. Based on our calculation, we propose that star formation in 
DLAs can trigger further star formation in dense clumps inside the 
DLA. 

When this cooling instability is triggered by the passing I- 
front, as the cloud collapses its H2 column density can temporarily 
reach very high values, up to 10 22 cm -2 for clouds with r ~ 5 pc. 
However, this phase is quite short-lived and the cloud rebounds and 
re-expands moderately within less than 1 Myr. Due to this short 
timescale and the very small cross-section of the clouds, we do not 
expect that such phase would be easily observable. On the other 
hand, smaller clouds, with initial radii < 4 pc, also collapse but 
never reach H2 column densities above ~ 10 18 cm~ 2 and most of 
the hydrogen is not in molecular form. 

Here we theoretically estimate the star formati on rate of a 
DLA b y using the cloud data from the simulation by Hirashita et al. 
120031) shown in Figure □ There are Na = 1092 clouds of 
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r > 4 pc per each kpc 2 in this simulation, corresponding to 
mean spacing between the clouds of x = 10 3 /N^ 2 — 30.3 
pc and mean volume per cloud of V = x 3 = 2.8 x 10 4 pc 3 . 
The total volume of the disk simulated in iHirashita et alJ 1 2003 ) 
is Vtot = 0.785 kpc 2 x 100 pc = 7.9 x 10 7 pc 3 , correspond- 
ing to total number of clouds iVtot = Vtot /V ~ 2800. Assum- 
ing a 4 pc cloud, and that 85% of the gas initially in the clump 
is converted into stars on a time-scale of 3 Myr, as our simula- 
tions indicate, this corresponds to a star formation rate per cloud 
SFR c i = 0.85M c i oud /3Myr « 70 M /Myr, and total SFR of 
SFR = SFR cloud N to t ~ 0.25 M yr" 1 kpc" 2 . This star 
formation activity produces a local UV field. Hirashita & Ferrara 
(2005) related the surface density of SFR, Esfr, and the UV field 
normalized to the Galactic value, \> as follows 

Esfr = 1.7 x 10~ 3 x M Q yr" 1 kpc" 2 . (12) 

Therefore, the SFR derived above produces a strong UV field cor- 
responding to x = 146. The SFR estimated above is significantly 
higher than the value derived fro m recent o bservation s of 0.01 — 
0.05 M B yr" 1 kpc" 2 <Wolfe et alJ l2003t IHirashita & Ferraral 
120051) . The assumption that all the clouds collapse at the same time 
could be unrealistic, and in reality, the star formation activity in 
DLAs could be regulated. On the other hand, the external UV back- 
ground, acts as a negative feedback in terms of star formation, as 
in our Case 3. Thus, we propose that the star formation activity in 
DLAs is quiescent (i.e. not starburst-like evolution) and strongly 
regulated by a combination of their internal radiative feedback and 
the external UVB field. 

In Cases 1 and 2, the velocity dispersion within the Fb-rich 
regions is a few km s" 1 or less , consistent with the values derived 
from H2 -detected components iLedoux et alT20 03). The evaporat- 
ing flow of ionized hydrogen produces a typical velocity width of 
~ 20 km s" 1 . The coexistence of molecular clouds and evapo- 
rating flows may be the reason for the observed com plex veloc ity 
structures in Fb-detected components (e.g. lPetitiean et all2002l) . 

At the interface of molecular clouds and evaporating flows, 
H2 suffers UV radiation and it is excited by UV pumping. This 
can be the reason for the observed high excitation temperature ob- 
served for H2 rotational levels J ^ 2 (Hirashita & Ferrara 2004 
and references therein). We predict that highly excited H2 has a 
higher velocity dispersion than H2 in low excitation levels, qualita- 
tively explaining some observations (e.g., an absorpti on system at 
2 = 2.595 in the sight line of Q 0405 -443: ILedoux et all2003l) . 



4 SUMMARY AND CONCLUSIONS 

We have examined if star formation activity in damped Lya clouds 
(DLAs) can be triggered in the presence of UV radiation. We have 
performed high-resolution hydrodynamics and radiative transfer 
simulations of the effect of radiative feedback from propagating 
ionization fronts on high-density clumps. We have considered two 
sources of UV radiation field to which high-redshift (z ~ 3) galax- 
ies could be exposed: one corresponding to the UV radiation origi- 
nating from stars within the DLA, itself, and the other correspond- 
ing to the UV background radiation. We have found that, for clouds 
with a radius r > 4 pc, the propagating I-fronts created by local 
stellar sources can trigger cooling instability and collapse of sig- 
nificant part, up to 85%, of the cloud, creating conditions for star 
formation in a timescale of a few Myr. The passage of the I-front 
also triggers collapse of smaller clumps (with r <C 4 pc), but in 



these cases the resulting cold and dense gas does not reach condi- 
tions conducive to star formation. Assuming that 85% of the gas 
initially in the clump is converted into stars in clouds with r > 4 
pc, and using the stat istics of the clouds in the simulation of an 
entire galactic disk by IHirashita et alJ 1200 3l) . we obtain a surface 
density of star formation rate of ~ 0.25 Mq yr kpc" 2 . This is sig- 
nificantly higher than the value derived from recent observations. 
On the other hand, the background UV radiation which has harder 
spectrum fails to trigger cooling and collapse. Instead, the hard pho- 
tons which have long mean-free-path heat the dense clumps, which 
as a result expand and essentially dissolve in the ambient medium. 
Therefore, the star formation activity in DLAs is strongly regulated 
by the radiative feedback, both from the external UV background 
and internal stellar sources and we predict quiescent evolution of 
DLAs (not starburst-like evolution). We believe these conclusions 
are fairly generic and are applicable also to molecular clouds in 
other types of galaxies, as long as the environmental conditions are 
similar, in particular in terms of significant presence of metals and 
dust. 
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